Quantum dynamics in single spin measurement 
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We study the quantum dynamics of a model for the single-spin measurement in magnetic- 
resonance force microscopy. We consider an oscillating driven cantilever coupled with the magnetic 
moment of the sample. Then, the cantilever is damped through an external bath and its readout is 
provided by a radiation field. Conditions for reliable measurements will be discussed. 
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I. INTRODUCTION 

Quantum state manipulation has become a fascinating 
perspective of modern physics. With concepts developed 
in atomic and molecular physics the field has been fur- 
ther stimulated by the possibility of quantum computa- 
tion. For this purpose a number of individual two-state 
quantum systems (qubits) should be manipulated in a 
controlled way. Several physical realizations of qubits 
have been considered in recent years, but nano-electronic 
devices appear particularly promising because they can 
be embedded in circuit and scaled up to large number of 
qubits. In addition to manipulation, single qubit mea- 
surement is an important and delicate question. For in- 
stance, a force resolution of attonewtons is required to 
detect qubits represented by single spins pj. Nowadays, 
magnetic resonance force microscopy (MRFM) is striving 
for this ultimate goal 2] . The use of cantilevers to detect 
small forces is at the heart of force microscopy in all its 
various forms [3(. Therefore, it is one of the most promis- 
ing techniques to also achieve single spin detection. 

Usually, small magnetic samples of unpaired electron 
spins are detected using a MRFM technique based on 
cyclic adiabatic inversion (CAI) 0, and the quan- 
tum dynamics of measurement for this method has been 
recently discussed in |3, |6|]. However, a new technique 
called "oscillating-cantilever-driven adiabatic reversals" 
(OSCAR) has been proposed and implemented in 0, 
where a sensitivity of about one hundred spins has been 
demonstrated. Very recently, this detection protocol has 
been improved and a sensitivity equivalent to a single 
electron spin has been finally achieved In the OSCAR 
protocol, the driven cantilever causes the adiabatic in- 
version of the paramagnetic moment of the sample. The 
backaction of this moment causes a frequency shift of the 



cantilever, which is then detected by a fiber-optic inter- 
ferometer. A classical description of the spin-cantilever 
dynamics in the OSCAR technique has been given in 
Ref. . The quantum dynamics of the spin-cantilever 
system has been instead studied in Ref. 0, where the 
master equation including the decohering effect of a high- 
temperature environment has been solved numerically. 
The present paper investigates the quantum dynamics 
of the whole system including not only the spin and the 
cantilever, but also the effective measurement device, i.e., 
the fiber-optic interferometer, which is here represented 
by a quantized radiation mode. The spin state manifests 
itself in the probability distribution of the phase of the 
field and our treatment allows us to take into account 
also the backaction of the measurement meter (the field 
mode) on the spin-cantilever system. Differently from 
jTo| . our approach is based on the exact analytical so- 
lution of the quantum dynamics which is valid in the 
adiabatic limit and for not too large spin-cantilever cou- 
pling. The paper is organized as follows. In Section [H] 
we introduce the model, then in Section lull we solve the 
quantum dynamics. In Section llVI we present the main 
results. Section IV is for concluding remarks. 
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II. THE MODEL 



We consider a model where a ferromagnetic particle is 
mounted on the cantilever tip. A permanent magnetic 
field B z points in the z direction, while a rf field of mag- 
nitude B\ rotates in the x — y plane in resonance with the 
spin precession around z axis at frequency ui s . A radia- 
tion field is then used to monitor the cantilever's position. 
We describe the cantilever interacting with a single elec- 
tron spin-1/2 as a harmonic oscillator of frequency u> c , 
and the radiation field as a single electromagnetic mode 
of frequency u) r . Then, the dimensionless Hamiltonian 
for such a system (scaled by the factor Twj c ) can be writ- 
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ten as 



H = a f a + u' r tfb + u' s S z 

+ e [S x cos (u> s t) + S y sin (u> s t)] 
- V2ri(a + a + )5 z - K,tfb(a + a f ) . 

Here a is the ladder operator of the cantilever 



2nu) r \ k r 



(1) 



(2) 



where P and Z are the momentum and z-coordinate op- 
erators of the cantilever, and k c is the cantilever spring 
constant. Moreover, b, s are the ladder operators of the 
radiation mode with [b, b'] = 1, u>' r = ui r /oj c , lo' s = ui s /lj c , 
and Si (i = x,y,z) are the usual spin operators obey- 
ing the commutation rule [S x ,S y ] = 2iS z and its cyclic 
permutations. The first row of Eq. 0) is the free Hamil- 
tonian of the system, the second row is the driving term 
of amplitude e given by 0, ITo| 



(3) 



where 7 e is the electron gyromagnetic ratio. Finally, the 
third row contains the spin-cantilever interaction (due to 
the spatiall y v arying magnetic field B z ) with coupling 
constant |^.IToj 



? = — 

2 V uj c k c 



dB z 



dz 



(4) 



and the optomechanical interaction (due to the radiation 
pressure force) with coupling constant k, which, for a 
cavity interferometer with length L is given by 



T 



2,k c 



(5) 



The Hamiltonian Q in interaction picture with respect 
to H = u>' r b^b + uj' s S z reads 

H = a + a + eS x - y/2r)(a + a f )5 2 - K,tfb(a + a f ) , (6) 

which is the same model Hamiltonian considered in |10| , 
supplemented with the optomechanical interaction with 
the readout radiation mode b. It is then convenient to 
apply the it/2 rotation around the y axis in the spin space 
R — exp [1^(^/4)] exchanging S x and S z , to obtain the 
transformed Hamiltonian 

H = RHR) 

= a) a + eS z - (\/2r,S x + nrfb^j (a + a f ) . (7) 

Apart from the additional optomechanical interaction 
term, Eq. Q is equivalent to the Hamiltonian of a two- 
level atom with energy separation 2e interacting with a 
quantized radiation mode with unit oscillation frequency 



represented by the cantilever z- motion The equiv- 
alent dipole-interaction term in Eq. J7J contains also 
counter-rotating terms, which are usually neglected in 
quantum optical situations (rotating wave approximation 
(RWA)) [l2|. This approximation is however not justi- 
fied in the present case. In fact, in typical MRFM sit- 
uations, the spin is much faster than the cantilever and 
adiabatically follows the cantilever motion, i.e., £ » 1. 
To be more precise, the OSCAR technique is usually op- 
erated under the adiabatic condition e 2 ^> rja m , where 
a m is the mean amplitude of the cantilever oscillations, 
and under the full reversal condition e <C rja m 0, • 
However, Ref. |lCj shows that the best signal sensitivity 
is achieved for partial adiabatic reversals, i.e. e ~ T]Ct m . 
Along this line, we consider the limit e ^> 1 + -q together 
with k<1, which allow us to treat the interaction part 
of Eq. (0) as a small perturbation with respect to the free 
Hamiltonian a^a + eS z . Thus, keeping the lowest nonzero 
order terms (first order in k and second order in rj), we 
end up with the effective Hamiltonian 



H, 



eff 



- „t 



a^a + tS z + x a >o>S z — nb^b{a + a'), 



(8) 



where the effective spin-cantilever interaction \ is given 
by 



X = 



16r] 2 e 
4e 2 - 1 



(9) 



The third term of Eq.lJBJ, describes a dispersive inter- 
action in which the spin and the cantilever energies 
are independently conserved and one induces a state- 
dependent phase shift on the other. In particular, the 
cantilever undergoes a frequency shift ±x depending 
upon the value of S z , which is just the frequency shift 
which is detected in the OSCAR technique by means of 
a phase sensitive measurement of the radiation mode. 
Therefore we expect that by measuring the probability 
distribution of the field phase one effectively measures 
the spin component S z , which corresponds to S x in the 
initial, non-rotated frame. This is consistent with the 
classical description of MRFM 9] , in which the cantilever 
always measures the spin component along the effective 
magnetic field which, in the strong adiabatic limit of very 
large e considered here, is essentially aligned along the x 
direction. 



III. QUANTUM DYNAMICS 

To describe a realistic situation we consider the can- 
tilever plunged in a thermal bath at equilibrium tempera- 
ture T and subject to a viscous force with damping rate T 
(with T/ui c = 7). This situation corresponds to an ohmic 
environment , and since for typical MRFM cantilevers 
lo c is of order of some kHz, it is N = kgT /Tiuj c -l/2> 1 
(ks is the Boltzmann's constant) even at cryogenic tem- 
peratures, so that environmental effects are satisfactorily 



3 



described in terms of the Caldeira-Leggett master equa- 
tion [LH , which is valid in the high temperature limit and 
is described by the superoperator 



Cp = ~-[a + a\{a-a\p}] 



7 ksT 

2 TlLOr 



\a + a\ [a + 



a),p\ 



(10) 



where p is the density matrix of the rotated system. 
Therefore, the complete system dynamics will be gov- 
erned by the following master equation 



dt P 



H e f f ,p 



Cp, 



(11) 



where, due to the adoption of the dimensionless Hamil- 
tonian of Eq. (JSJ), time t is also dimensionless, i.e., scaled 
by uj c - 

In order to solve the master equation l|ll|) we introduce 
the eigenstates |e), \g) of S z with eigenvalues +1, —1 
respectively, and we split the total density operator as 



p {ee) ®\e){e\+p^ ®\g){g\ 

P {e9) ®\e){g\ + P {9e) ®\g){e\ 



(12) 



With p ( S^ = [p( e 9)] f . 

The terms p^' = e,g) in Eq. lfT2^l can be ex- 
panded on the set of Fock states external product op- 
erators |n) r (m| (n,m £ N) for the radiation mode and 
on the continuous set of displacement operators D C (X) = 
cxp{Aa^ — A*a} (A <E C) for the cantilever mode, as follow 

m 

? w) (t) 



n,m— 



d 2 x 



— e"l A l / 2 eW)(A,t)£>t( A )®|n) r (m|, 

(13) 

where Qnm (t) represent normally ordered characteristic 
functions |l5l Il6j | for the cantilever mode for given spin 
and radiation states (they are explicitly given in the Ap- 
pendix [SJ . 

From Eq. (|12f> one gets the non-rotated solution in the 
original frame 



p(t) 


= R^p{t)R 






(14) 


1 

~ 2 


p^(t)+p^{t) 


-p^\t) 


-p 


®|e)(e| 


1 

+ 2 


~p( ee \t)+p( gg) (t)~ 


¥p {eg) {t)- 


¥p {9e \t) 


®\g)(g\ 


1 

+ 2 


p (ee) {t)~p {99) (t)- 


¥p {eg \t)- 




®\e){g\ 


1 

+ o 


~p {ee \t)-p {9a \t)- 




¥ P (9e) {t) 


®\g){e\. 



As initial condition we consider a coherent state \/3) r for 
the radiation field, while for the cantilever we consider 
a thermal state at temperature T, shifted by a complex 
amplitude a, where Re{a} gives the initial mean value 
of the position and Im{a} the initial mean value of the 
momentum. 



IV. QUANTUM MEASUREMENT 

As discussed above, the cantilever's position carries in- 
formation about the spin state and, in turn, it is read by 
the radiation field. Due to the optomechanical coupling, 
the cantilever affects the phase of the field and we have 
therefore to look at the probability distribution of the lat- 
ter. The canonical probability operator measure (POM) 
for phase measurement is [l7| 



11(6) 



27T 



with 



E 



e in6 \n) r . (15) 



Hence, the phase probability distribution results 
P{6 lT ) = Tt{IL(9)p(t)} = P ee (9,T)+P gg (9,T) 



(16) 



where we have introduced the new time r = which, 
due to the scaling chosen above, is a dimensionless time 
now scaled by the cantilever damping rate F, and 

P jj (e,r) = Tv{ll(0)p m (r)} , (17) 
with j — e,g. From Eq. (|13|1 we get 

1 oo 

e-^-™) 9 eto)(A = 0,r). (18) 



n,m— 



Let us now determine the conditions under which the 
apparatus performs reliable measurements of the spin. 
As discussed above, the apparatus is expected to mea- 
sure S x (in the original non-rotated frame) and therefore 
the easiest case is when the initial state of the spin is 
an eigenstate of S x (in the non-rotated frame). We con- 
sider, for example, the eigenstate with S x = —1, which 
correspond to state \g) in the rotated frame and therefore 
implying 



e(*»(A,r = 0)=e- 
x exp < a* A — a* A - 



V n!m! 

|A| 2 W (19) 



Twj, 



As a consequence, S z = — 1 (in the rotated frame) is 
a constant of motion and one has a cantilever frequency 
shift x (see Eq. (|SJ)) which manifests itself as a shift of the 
phase distribution (see Fig. [TJ. In fact, the phase proba- 
bility distribution is initially peaked around 8 = 0, due to 
the real amplitude (j3 = 3) of the initial coherent state of 
the radiation mode. Then the peak moves, toward larger 
values of 9, and faster than that with x = 0. If we had 
started from state |e) (S z = 1 in the rotated frame), one 
would have a smaller phase shift of the peak. This could 
be surprising as one would expect a greater phase shift for 
frequency (1 + x) an( i a smaller one for (1 — x) ■ However, 
if one just considers the unitary evolution arising from 
Eq.©, one reduces to the model studied in Refs.[Tl||. 
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FIG. 1: Phase distribution for an initial state of the spin 
(\g) — \e))/y/2 (in the non-rotated frame) at times r = 
(dotted line), r = 8 x 10 4 (dashed and solid line). The dashed 
line refers to \ — while the solid line to \ — 0.5. The values 
of other parameters are 7 = 10 -4 , /t/7 = 0.08, a = 4i, (3 = 3, 
and N = 10 2 a), N = 10 4 b). 



FIG. 2: Phase distribution for an initial state of the spin 
\e) (in the non-rotated frame) at times r = (dotted line), 
r = 8 x 10 4 (dashed and solid line). The dashed line refers 
to x = while the solid line to x = 0.5. The values of other 
parameters are 7 = 10 -4 , ft/7 = 0.08, a = 4i, ft = 3, and 
JV = 10 2 a), JV = 10 4 b). 



where it is possible to argue that the induced phase shift 
is inversely proportional to the cantilever's frequency. 

In Fig. n we have considered two different values of 
temperature, N = 10 2 (a) and N = 10 4 (b), which for 
typical cantilever frequencies uj c /2ir ~ 10 4 , correspond to 
T ~ 60 /iK and T ~ 6 mK respectively. One sees that 
already at mK temperatures the peak in the phase distri- 
bution corresponding to the state of the single spin starts 
to become too wide and flat to be clearly distinguished. 

The conditions which the various parameters have to 
satisfy in order to have an unambiguous discrimination 
of the two spin eigenstates can be better established if we 
consider a different initial condition for the spin, i.e. an 
initial superposition of the two eigenstates of S x (in the 
non-rotated frame). Fig. prefers to the initial condition 
S z = 1 in the initial frame, corresponding to the state 
(\g) + |e))/v2 in the rotated frame. This implies 



eW)(A,r = 0) = e- 



\Jn\m\ 



xexp<ja*A-a*A-^|A| 2 
huj c 



j = e, g.(20) 



Due to the superposition, the measurement performed 
by the MRFM is reliable when we have the simultaneous 
presence of two distinct peaks in the phase distribution, 
associated with the two different phase shifts induced by 



the two spin components. This means that the measure- 
ment time has not to be too small, in order to allow the 
two peaks to become distinguishable. This implies the 
condition \t = x( T /l) > 1- Then, reliable measurement 
can be performed until the damping affects the dynam- 
ics, as it is confirmed by Fig. where one still has two 
distinct peaks corresponding to the two spin states when 
7T = 8. In Fig. the state of the single spin is well 
detected because we consider the low temperature con- 
dition N = 10 2 , i.e., T ~ 60 /LtK (see above). At higher- 
temperatures, the effect of the thermal environment is 
instead more destructive, as it can be seen from Fig. ^jp, 
which corresponds to N = 10 4 (T ~ 6 mK), where the 
two peaks become too broad and are no more distinguish- 
able. 

The other parameter values in Figs. 1 and 2 are x = 0.5 
and 7 = 10~ 4 . This value of 7 is realistic, because 7 -1 
corresponds to the cantilever's quality factor Q, which is 
even higher in recent experiments (see Ref . 8] ) , while the 
chosen value of x is still far from present experimental 
values. In fact, from Eq. © one has x — 47y 2 /e and using 
Eqs. © and Q and the experimental values of Ref. 0, 
one gets e ~ 800, 77 — 0.04, so that \ ~ 1°~ 5 . Unfor- 
tunately, for such a small value of Xi the radiation field 
phase shift due to a single spin is too small to be detected 
and therefore larger values of Xi which means larger val- 
ues of the coupling constant 77, are needed to allow a 
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P(6) 



1 . 2 



. 8 
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the presence of a high-temperature ohmic environment 
acting on the cantilever. In this limit, the probability 
distribution of the phase of the optical readout mode 
gives a clear signature of the spin component along the 
effective magnetic field provided that the effective spin- 
cantilever coupling coefficient \ of Eq. 10 is sufficiently 
large and the optomechanical coupling k is either not too 
large (in order to limit the backaction) or not too small 
(in order to extract information from the system). 
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FIG. 3: Phase distribution for an initial state of the spin 
e) (in the non-rotated frame) for different values of k at time 
r = 8 x 10 4 . It is k/7 = 0.04 dotted line, K/7 = 0.08 solid line 
and k/~/ = 0.12 dashed line. The values of other parameters 
are 7 = 10" 4 , x = 0.5, a = 4i, f3 = 3, N = 10 2 . 
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APPENDIX A 



reliable detection of the spin state within a reasonable 
detection time. Also a larger number of photons in the 
radiation mode would improve the resolution because in 
such a case the peaks in phase distributions would be 
narrower. 

Finally, another important parameter is the optome- 
chanical coupling k between cantilever and the radiation 
readout mode. Similarly to %, this parameter has to be 
sufficiently large, because otherwise the phase measure- 
ment has not enough sensitivity (see the dotted line in 
Fig. |3 where the two peaks are not well resolved) . On 
the other hand, the readout mode is a quantum system 
and its backaction can blur the measurement result. This 
happens if the coupling k is too large, as it is shown by 
the dashed line of Fig. [3] In this case the backaction is 
responsible for a phase diffusion of the peaks which are 
washed out and no more distinguishable. By inspection 
from Eq.JHJ) the condition k\(3\ 2 <C \a\ 2 should guaran- 
tee a negligible optomechanical backaction. However, it 
must be modified to take into account the elapsed time 
whenever the measurement is performed after many cy- 



cles, which naively implies kt\ 
the values in Fig.^J). 



<C I a 1 2 (confirmed by 



Since only 6™ and Qnm enter in Eq. (|TtJ)l . we need 
to solve Eq. (fTTf) for and p^g) only. 



1. Solution for p (e 



Let us define pnm = r(«|p^ ee ^| 
be rewritten in the form 



then Eq. l(TT|l can 



dee) 



dr' 



■ 1 + X 

7 



" "i Cum 



+ i-n(a + at)pj$ - i-mp^ (a + a f ) 
7 7 
1 



+ ijV (2otp - aatpM - p^aa^) 
- \M (2atpM a t _ ( a t) 2 - ?H ( a t) : 



-M * [ 2ap^e) a _ a 2~(ee) _ ~(ee) a 2 



(Al) 



CONCLUSIONS 



with 



We have presented a model for the description of the 
quantum dynamics of the measurement of a single spin 
by a MRFM operating with the OSCAR technique 0- 
Our treatment includes also the phase-sensitive readout 
of the cantilever frequency-shift provided by a quantized 
radiation mode describing the measurement performed 
by an optical interferometer. We have solved analytically 
the dynamics in the strong adiabatic limit of a very fast 
spin, in which the spin reversal cannot be complete, in 



N : 



kuT I 



Twj c 



M 



k B T 



(A2) 



and r = -ft. Eq. i|AI|) is a valid master equation, preserv- 
ing t he p ositivity of the density matrix for 7V(7V + 1) > 
|M| 2 [T{| which occurs in the high temperature limit, i.e. 
ksT ^> hu> c . 

Since Qnm represent normally ordered characteristic 
functions, we can transform Eq. IjAIfl into a partial dif- 
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ferential equation [T^, that is 
.1 + X 1 



7 

1 + X 

7 



9A 



<9A* 



k d k d 

-i—vn — m)— — i—[n — m)—— 
7^ ' dX Y ' dX* 

-i—mX — innX* 

7 

1.9 1,. 



2 X 'dX* 
N 



2 X *dX 



N 



x 2 - ^ (x*f - mxf eg 



(A3) 

We then use a Gaussian trial function, which is justified 
by the choice of an initial Gaussian state for the cantilever 
(a shifted thermal state) and by the fact that Eq. (|A3() 
preserves Gaussian states, 

e^(r) = e^{A(T)X + B(r)X* +C(r)\X\ 2 



(A*r 



+ F(r) . (A4) 



where we have set a; = (l+x)/7- Accordingly to Eq. 119|l . 
initial conditions now read 



A(0) 


* 

= a 


5(0) 


= — a 


C(0) 


--( 


D(0) 


= o, 


E(0) 


= o, 


F(0) 


= In 



Notice that 



/3" (/?* 



,-|/3|- 



V 'n\m\ 



(All) 
(A12) 

(A13) 

(A14) 
(A15) 

(A16) 



As a consequence, Eq. (|A3|I yields the following set of 
linear differential equations 



A(r) 



B(r) 



1 



1 



iu --)A(r)--B(T) 



-i-(n - m ) (CM - DM) - i-m , (A5) 

7 7 



1 



AM 



bm 





■ K i 

+i — [n 


- m) (CM - EM) - i-n 


(A6) 




7 


7 


cm 


= -CM 


-1(E(t)+D{t))-N, 


(A7) 


dm 


= -CM 


+ (2iw- I) DM ~N, 


(A8) 


em 


= -CM 


- (2iuj+l)EM ~N, 


(A9) 


fm 


.K 

= t — (n — 

7 


m) (AM BM) , 


(A10) 



9(-)(A = 0,r) = Tr ( 



{,<->} 



Solving Eqs. l|A"5)) - l|A10|) we get 



(A17) 



FM = F(0) + 



F, 



i + n 



i(n/"/)(n — m) 
4j 2 - 16w 4 



3 -(1+0)t _ 1 



-2 



Fi 



i + n 

Fi 



i - n 



e -(l+H)r/2 _ x 



- 2 



1 - ft 



F 5 [e 



- 1 



e _(l_n)r/2 _ 1 
I +^6T 



(A18) 
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where il 2 = 1 — Aui 2 and 
F 1 = 



F 5 
F e 



— (n — m) 

7 



(J 2 LU N 

AiN + l6iNoj 2 T 8i— (1 + 3N) ± 32i— (1 + N) ± Ai — 



H — (n + m) 

7 



-2w + 8u> ± 2— =F 8— 



+Re{a} 



4w 2 ± 16w 4 + 4— ± 16— 



Im{a} 



UJ 3 Ll> 51 



7 



■(n — m) 



n 1 

-2i + 4iuj 2 ±2i-T8i— 



8iu —(n — m) , 

7 

-(n-m) [4i + 8iN - 16iw 2 - 32iw 2 7V] + -(n + m) [4d - 16w 3 ] 



r 



(A19) 

(A20) 
(A21) 
(A22) 



2. Solution for p (99) 

The quantities 8^? (A = 0, r) are given by Eas. (|A17l) - 
22|) with the replacement x ~* ~X (hence (1 — x)/l ~ ¥ 
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